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We study a one dimensional gas of repulsively interacting ultracold bosons trapped in a double¬ 
well potential as the atom-atom interactions are tuned from zero to infinity. We concentrate on the 
properties of the excited states which evolve from the so-called NOON states to the NOON Tonks- 
Girardeau states. The relation between the latter and the Bose-Fermi mapping limit is explored. 
We state under which conditions NOON Tonks-Girardeau states, which are not predicted by the 
Bose-Fermi mapping, will appear in the spectrum. 


The recent ground-breaking experiments on the trap¬ 
ping of a few fermions or bosons have stir up the the¬ 
oretical interest in small ultracold quantum gases [H-Q. 
The atoms in these experiments can be often regarded 
as being effectively trapped in one dimension and inter¬ 
acting through contact interactions. This is a system 
with many examples of exact solvability (see, e.g., 0“ 
i). A step forward in these experiments is to consider 
other geometries, a goal which is already accomplished 
for double-well potentials 0 - 


Experiments have already reached the strongly inter¬ 
acting ID gas of bosons 13, 1^, exploring the mapping 
of repulsive bosons into a system of ideal fermions to 
form the so-called Tonks-Girardeau (TG) gas 15, 0. 
Simultaneo usly , double-well experiments with ultracold 
bosons 171, ll8| should sooner than later allow one to pro¬ 
duce macroscopic superpositions 0 |. So, combining the 
two may allow one to produce macroscopic superposi¬ 
tions of strongly interacting TG gases. In this letter we 
settle under which conditions the superpositions in the 
Fock regime, so called NOON (cat-like) states 2nl-22|. 
evolve into superpositions of TG gases, termed NOON 
TG, as the interactions are tuned from small to infinity. 


Model Hamiltonian. Let us consider N atoms 
trapped in a one-dimensional double-well potential, 
described by the first-quantized Hamiltonian H = 


E;=i [-(1/2)52,+Fdw(^,) 

tion is taken as a contact, 14„ 


+ Fnt, where the interac- 

int = - Zj) and 

where h = m = 1. To model the double-well poten¬ 
tial we use a simplified Duffing potential, Vuw[z) = 
Vb (—8(z/Vb)2 -b \&{z/VqY + l)- This Duffing potential 
has two minima equal to zero at z = dc.VQl2 and a local 
maximum equal to Vq at z = 0. Vq is thus the bar¬ 
rier height and the distance between the wells. In our 
study, we take Vq € [Iq™™, with Vb™™ a small 

number and large enough to get a degeneracy of 

the order of I0“^° between the first two single-particle 
eigenstates (that is ~ 12). By increasing Vq, the 

confining potential goes from a quartic, for Vq —>■ 0, 
into increasingly deeper double-well potentials. Let 
be the n-th eigenstate of the single-particle Hamiltonian 



FIG. 1. (a) The first ten eigenvalues of the single-particle 
Hamiltonian as a function of Vb are depicted. The dashed line 
represents the condition E = Vo- (b)-(d) One-body density 
matrix for NOON Tonks-Girardeau states of A = 2, 3 and 4 
obtained with Vb = 6,10 and 12, respectively. 


Hsp = (1/2)52, -b 14)w(z) with energy e„. As Vq is 
increased, quasi-degeneracies between the single-particle 
energies are introduced, see Fig. [Ua). These occur for 
values of the corresponding quasi-degenerate pair of en¬ 
ergies smaller than Vq. 

Spinless fermions. For N spinless fermionic atoms, 
the ground and all excited states tVl of H are Slater de¬ 
terminants built up from N single-particle eigenstates />„. 
This wave functions fulfil the Pauli exclusion principle, 
i.e. are zero whenever Zj = Zj, and therefore do not feel 
the contact interaction potential. Also, these solutions 
have a sign change at these points, i.e. for fermions, the 
wavefunction has to change sign when an atom of coor¬ 
dinate Zi is interchanged with one at zj. The energies 
of the many-body eigenstates are the summation of 
the energies of the single-particle states used to build the 
corresponding Slater determinant. 

For N = 2 fermions, which will be used as 

a limiting case in the forthcoming discussions, the 
ground state is 'I'g(zi,Z 2 ) = l/-\/2Slater(())o,^!)i) = 
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FIG. 2. Energy spectra for A'" = 2 (a) and = 3 (b) with 
Vb = 5 and Vo = 10, respectively. The blue dashed lines in 
both panels represent the energies expected in the Fermion- 
ized limit, g ^ oo (we explicitly label a few of them to help to 
understand this limit). These energies have been calculated 
by extrapolating the results obtained with an increasing num¬ 
ber of modes, with a maximun number of modes of 30 and 24 
for N = 2 and N = 3, respectively. 


l/\/2{(j)o{zi)4ii{z2) - (/'i(zi)(/)o(-Z 2 )) and has energy £^0 = 
£0 + £1 for all Vq. For large Vb, the first two single¬ 
particle states are quasi-degenerate, ei = £0 + £ 01 , with 
£01 -C £ 0 . The next four excited states are 
Slater((()o, <(' 2 ), ^2 Slater((()i, (/) 2 ), d/l ~ Slater((()o, 
and ~ Slater((()i, i/fa). For Vq > 4.5, the second and 
third single-particle states are also quasi-degenerate (see 
Fig. m, and therefore their energies can be written as 
£3 = £2 +£ 23 , with £23 “C £ 2 - Thus, for Vq > 4.5 the first 
four excited states are quasidegenerate and have energies 
£1 = £0 + £ 2 , E2 = £0 + £2 + £ 01 : = £0 + £2 + £23: and 

£4 = £0 +£2 +£01 +£ 23 - For two atoms, the next excited 
state has a finite energy gap with this manifold (see gap 
between £3 and £4 for large Vq in Fig.[T]). 

For sufficiently large Vb, such that a sizeable amount 
of pairs of single-particle states are quasidegenerate, we 
can define the following single-particle wavefunctions. 


i>n — —7={4'2n + (t>2n+l) 

W = 2 

i'n = - 4>2n+l) 


,n = 0, 1, 


( 1 ) 


Here j = L{R) stands for left (right). These functions 
are mostly localized either in the left or right well when 
the pair of delocalized functions used to construct 
them are quasi-degenerate. For large interaction ener¬ 
gies, the small degeneracy breaking between the states 
will in some cases become irrelevant, and a much simpler 
picture is obtained in the localized basis. 

Interacting bosons in a double well. Using the 
single-particle eigenfunctions to build the Fock basis, the 
many-body Hamiltonian for bosons reads 


M M 

■H = ^£iajai-I-I ^ Ik,i,m,na\a\aman, (2) 

i—1 k,l,m.n 



Z Z 


FIG. 3. One-body density matrices of two of the four excited 
states for two bosons in a double-well for Vo = 5. The upper 
(lower) row is computed with g = 0.5 {g — 20). One-body 
density matrices for a NOON state, panel (a), and its evolved 
state, panel (c), and a non-interacting state, panel (b) and its 
evolved state, panel (d). 


with = J dz4>k{z)(j)i{z)<prn{z)4'n{z), whcrC 

[d\^di] = Sik and [d[,,aj] = [afc,ai] = 0. A Fock vector is 
written as |no, ni,... um) = A/’(aJ)”“ ..., (aj^^)”" |vac), 
where |vac) is the vacuum, rii is the number of atoms in 
mode i and = \/no!... um-- 

Beyond the bimodal Fock regime. For small in¬ 
teractions the system remains bimodal, allowing to ex¬ 
plore with a two-mode model from the Josephson to the 
Fock regimes 11, 1^. As interactions are increased, the 

a item approaches the strongly correlated regime 23- 
. Let us explore in detail this transition, using as 
driving examples the N = 2 and N = 3 cases. 

Assuming a barrier high enough to have two quasi¬ 
degenerate pairs in the lower part of the single-particle 
spectrum, the interaction energy will make the energy 
gaps irrelevant. In this case, the simplest picture 
is provided using the localized single-particle ba¬ 
sis defined in Eq. (HI)- Using this basis, the Fock 
vectors can be written as, jug , n^, nf, nf-,... )io = 
Ar(ao ^)"o (a^ + )"? ... (a^^/ 2 )"“'''(«M/ 2 )”"^'|vac) with 


N ^ = ,/n^!ng! 


T? ^ ^77 ^ ^ 

■ • "'M/2-'*M/2" 


For instance, in the N = 2 case, once the interac¬ 
tion is larger than the gap, the ground state is well ap¬ 
proximated by |1,1,0,... )io. This state has one atom 
in each well and is thus unaffected by the interac¬ 
tion. Its energy remains constant as g is increased, 
see Fig. [DJa). The first two excited states, quaside¬ 
generated with the ground state in the non-interacting 
case, are the NOON states, INOON 2 ) = (|0,2,0,... )io ± 
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FIG. 4. Comparison of two of the four excited numerical wave 
functions for g = 20 (upper row) and the analytic ones Vlti 
(e,g) and ^'3 (f,h) obtained through the Bose-Fermi mapping 
(lower row). Panels (a), (b), (e) and (f) are obtained with 
Vo = 4, while panels (c), (d), (g) and (h) correspond to Vb = 5. 
N = 2. 

|2, 0, 0,..., )io)/-\/2, which contain a delocalized pair of 
interacting atoms (two atoms in the same well). Their 
energy grows linearly with g for small g, see Fig. [2](a). 
The next excited states involve more than two modes 
and are clearly gapped in absence of interactions [ 2 ^, . 

They read, (|0,1,1,0,... )io± |1,0, 0,1,0,... )io)/^2, and 
(|0,l,0,l,0,...)io± |1,0,1,0 ,...)io)/a/2. As before, the 
first two are mostly non-interacting and thus their ener¬ 
gies remain constant as g is increased. The latter, how¬ 
ever, have one pair of interacting atoms, their energy 
grows linearly with similar slope as the INOON 2 ) states, 
see Fig. [2Ka). 

For N > 2 there are no noninteracting states. For 
N = 3 the lowest energy manifold at small interactions 
is (|1, 2, 0 ,... )io ± |2,1, 0 ,..., )io) and the NOON states, 
INOON 3 ) = (|3, 0 , 0 , ... )io ± | 0 , 3, 0 , ..., )io). The pair of 
atoms on the same well in the first ones are responsible 
for the linear increase in the energy as a function of g, 
with same slope as for TV = 2 discussed above. In con¬ 
trast, the NOON states have 3 interacting pairs, with a 
correspondingly larger slope, see Fig. [2Kb). 

In the general case, the lowest energy manifold at small 
g is |iV,0,0,...)io± |0,Af,0,...)io, |1,V- l,0...)io± 
|7V- 1,1,0...) lo, etc. Our interest at this point is to 
disentangle the fate of these states as the interaction is 
varied from the Fock regime into the strongly interact¬ 
ing regime. In particular, what happens with the states 
in which more than one atom populates each well in the 
Fock regime. The driving intuition is twofold. On one 
hand, we know that for bosons in a single well, the sys¬ 
tem evolves into a TG gas, i.e. the NOON would directly 
evolve into a NOON TG gas. On the other hand, the 
Bose-Fermi mapping theorem can be directly applied to 
the system, providing exact solutions in the infinitely in¬ 
teracting case, which will be shown to be in contradiction 
with the first one. 

Prom NOON to the NOON Tonks-Girardeau. 

The N = 2 and N = 3 cases are again very illustra¬ 


tive. We note that energetically the fate of the NOON 
states in iV = 2 is very similar to what happens to the 
(|1,2,0,... )io ± |2,1,0,... )io)/-\/2 states in N = 3, see 
Fig. 121 i.e. the third particle in the N = 3 plays an 
spectator role as g is increased. The latter takes place 
also in the general N > 2 case, where none of the eigen¬ 
states is non-interacting and in which upon increasing g 
the system is expected to increase correlations to avoid 
the interaction. Thus, let us first disentangle the fate 
of the NOON states at N = 2. From Fig. |2Ka) we see 
that as g is further increased, their energies saturate to 
a constant value which actually approaches £0 +£ 2 - This 
is also essentially the energy of the states evolved from 
(|0,1,1,0,.. .)io ± 11,0,0,1,0,.. .)io)/-/2, which are non¬ 
interacting. Thus, in the strongly interacting regime, 
there are four quasi-degenerate states. 

The evolution with increasing g of their one-body den¬ 
sity matrices (OBDM) is very telling. The NOON states, 
evolve from the INOON 2 ) states at g = 0.5, see Fig.jSKa), 
to distributions with two peaks per well (or three and four 
peaks per well for N = 3 and N = A, respectively, see 
Fig. |TKc,d)), resembling the OBDM for two TG gases, 
Fig.lSIc). The OBDM for the non-interacting states sim¬ 
ilarly shows two-peaks per well, as would correspond to 
populating the single-particle states (f> 2 , 3 , which have one 
node in each well. 

Let us now examine the prediction from the Bose- 
Fermi mapping theorem. The rigorous Bose-Fermi map¬ 
ping theorem is established by noting that the bosonic 
problem in the ^ 00 is equivalent to the fermionic 
problem when imposing the boundary condition that 
the wave functions have to vanish if Zj = Zji. This 
boundary condition is obeyed by the fermions for free 
because it is imposed by the Pauli exclusion principle. 
Then, the bosonic wave functions can be obtained 
from the fermionic ones after symmetrizing them, 
with A = nj>j'Sgn(Zj — zj/), with sgn(z) 
the sign function. The analytic form of the wave func¬ 
tions for the case of fermions was discussed above. The 
Bose-Fermi mapping theorem applies equally regardless 
of the value of Vq. 

In Fig.lHwe compare the predictions of the Bose-Fermi 
mapping, to the actual numerical results, for N = 2 
and for Vq = A and Vq > 4.5. For Vb = 4, the Bose- 
Fermi mapping, panels (e,f), provides an accurate rep¬ 
resentation of the numerically obtained results, panels 
(a,b). Both the NOON and the non-interacting states, 
evolve at g = 20 to the states predicted by the Bose- 
Fermi mapping. Note that in this case, the barrier is not 
high enough to substantially localize the single-particle 
ground and first excited state, see Fig. |TJ None of the 
many-body states is however a NOON TG state. 

A surprising result appears however as we increase the 
barrier height, going from Vb = 4 to Vb = 5 in Fig. [H In 
this case, the prediction of the Bose-Fermi mapping, pan¬ 
els (g,h), and the numerical results, panels (c,d), clearly 
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FIG. 5. Detailed analysis of the g —>■ oo limit for Vo = 5.5 
(upper panel). Two of the first excited states for N = 2 and 
Vo, computed numerically, solid lines, compared to the results 
of the perturbation theory when the width of the perturbation 
is cr ~ 0.07 (see the text for further details), dashed lines. The 
small panels depict the composition of the dressed states, 4'i 
and (4'fc|4t5) (solid, red) and ('J'llT^) (dashed, blue). A 
similar qualitative picture is obtained for states 4'2 and 4'|. 


disagree. The numerical results show less left-right co¬ 
herence than the Bose-Fermi prediction. That is, there 
is essentially zero probability of finding two bosons in the 
same well, as seen in Fig. |4] (d), while in the Bose-Fermi 
mapping, the probability is clearly non-negligible in all 
four states. It is worth mentioning that these differences, 
which involve pair-correlations, are not reflected in the 
density profiles, i.e., the diagonal of the one-body den¬ 
sity matrices. Figs. [TKb,c,d) and[3](c,d). What causes 
such discrepancy? And, more important, what would 
an experiment find? Let us note, that in all cases, the 
left-right symmetry is fully respected, i.e. it is not a con¬ 
sequence of spurious numerical biases in the numerics. 

Finite-range effects. The Bose-Fermi mapping ap¬ 
plies for delta interaction potentials in the strict infinite 
interaction case. The Slater determinant ensures that 
the atoms do not interact, but, any finite size correction 
to the delta contact potential immediately implies a non¬ 
zero interaction between the two atoms. In addition, any 
numerical approach has an inherent minimal distance re¬ 
solved, in our case related to the maximun number of 
modes considered. Thus, for any number of modes, if g 
is sufficiently large, finite size corrections to the contact 
should show up. This energy stemming from the finite 
range of the atom-atom interaction, can be comparable 


to the splitting between the four excited states described 
above, resulting in a mixture of the four states. 

This hypothesis can be tested quantitatively. Let 
us consider a finite range perturbation on the four 
states predicted by th e Bo se-Fermi mapping, e.g. 
Vpertizj - Zji) = rjjV^TT^ exp[-{zj - 
and diagonalize it in the restricted space spanned by 
\]/b vj/b v]/b} Then we obtain the dressed states 
{vjib ilfb vjyb vjybj perturbation theory prediction 

can be compared to the numerical results by noting that 
the rj = 0 corresponds to the value of g for which the 
bose-fermi energies are reproduced for the number of 
modes used. As seen in Fig. [5] the model predicts a 
nice avoided crossing in agreement with the numerics. 
For r] < 0, thus removing repulsion from the g —>■ oo 
limit, we have the NOON TO as the lowest states of the 
manifold, — (^i + (see Fig. [S]) and ^2 — 

For 77 > 0 , the repulsive finite size correc¬ 
tion sets the non-interacting states — ('^'i ~ 
and ^2 — (^2 ~ as the lowest in the manifold. 

These results, checked for N = 2,3 and 4 are expected 
to remain valid for larger number of atoms. NOON 
states will evolve to NOON Tonks-Girardeau states pro¬ 
vided two conditions are met: 1) the potential barrier has 
to be large enough to ensure the existence of N quasi¬ 
degenerate doublets in the single-particle spectrum, and 
2) the actual range of the contact interaction is non-zero 
but finite, which is what actually happens in an exper¬ 
iment, such that in the infinitely interacting limit the 
residual interaction mixes the corresponding states pre¬ 
dicted by the Bose-Fermi mapping to produce NOON 
Tonks-Girardeau states. These findings, which require 
two-body correlations to be explicitly measured, will 
be of relevance in future experimental efforts to build 
strongly correlated states with inter-site spatial entangle¬ 
ment. As an outlook, we foresee that, for a general trap¬ 
ping potential with an associated single-particle eigen- 
spectra showing quasidegeneracies, the effect of the finite 
range of the contact interactions may involve discrepan¬ 
cies between the result predicted by the Bose-Fermi map¬ 
ping and the actual solution. 
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